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ABSTRACT 

We have performed the detailed numerical simulations on the propagation of 
the UHE protons in the energy range E = (10^^'^ — 10^^ °) eV in the relatively 
strong extra-galactic magnetic field with strength B = (10, 100) nG within about 
40 Mpc. In this case, the deflection angles of UHECRs become so large that 
the no counterparts problem is simply solved. As for the source distribution, we 
assumed that it is proportional to the number distribution of galaxies within 
the GZK sphere. We have found many clusters, which mean the small-scale 
anisotropy, in our simulations. It has been also shown that the observed energy 
spectrum is well reproduced in our models without any fine-tuned parameter. 
We have used the correlation value in order to investigate statistically the 
similarity between the distribution of arrival directions of UHECRs and that 
of galaxies. We have found that each correlation value for each parameter set 
begins to converge when the number of the detected events becomes O(IO^). 
Since the expected number counts by the experiment of the next generation such 
as TA, HiRes, Auger, and EUSO are thought to be the order of 10^, we will be 
able to determine the source distribution and values of the parameters in this 
study in the very near future. Compared with the AG AS A data, the significant 
anisotropy on the arrival directions of UHECRs are found in the analysis of first 
and second harmonics. This may originate from the incompleteness of the ORS 
database. This problem may be also solved if the source distribution is slightly 
changed. For example, this problem may be solved if we assume that UHECRs 
come from some of the galaxies such as AGNs and radio galaxies. 

Subject headings: cosmic rays — methods: numerical — ISM: magnetic fields — 
galaxies: general — large-scale structure of universe 
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1. INTRODUCTION 



The observed differential cosmic ray spectrum is remarlcably featureless and extends 
beyond 10^° eV ([lakeda et al. 1999| ). So far, the observed number of the Ultra High 
Energy Cosmic Rays (hereafter UHECRs) whose energies are above 10^° eV (this is the 
definition of UHECRs in this study) is only 25 ([Virmani et al. 2000|) . On the other hand, 



a new generation of the ground-based large aperture experiments such as Telescope Array 
(hereafter TA; P?elescope Array Project Design Report 2000|) , HiRes ([Wilkinson et al] 



|1999|) and South and North Auger QCapelle et al. 19^ ) is expected to detect about 1000 



UHECRs until 2010 ( |Zas 2001] ). Moreover, the mission proposed as EUSO is devoted to the 



exploration of UHECRs from satellites and is expected to detect about 1000 UHECRs per 
year ( [Bhattacharjee and Sigl 1998| ). We will be able to discuss statistically the feature of 
UHECRs, such as arrival direction, in the very near future. 

The mechanism to produce such highly energetic cosmic rays is still controversial. 
In most of the conventional acceleration scenarios, which are called bottom-up scenarios. 



effects of diffusive shock acceleration are taken into consideration (e.g., [Biermann 1995 



Halzen and Zas 1997| ; [Waxman 2000| ). So far, by using Hillas-plot ( pRllas 1984] ; [Selvon 



pOOO|) , gamma-ray bursts (GRBs) and/or active galactic nuclei (AGN) are considered as 
probable candidates. On the other hand, there is a large number of production models 
based on exotic particle physics scenarios (see [Bhattacharjee and Sigl 2000| and references 
therein). Sometimes effects of interaction or collapse of Topological Defects are taken into 
consideration to produce such massive exotic particles (e.g., [Bhattacharjee et al. 199^ ). 
These scenarios are called top-down scenarios. In this study, we mainly consider bottom-up 
scenarios although top-down scenarios are very attractive and important. This is because 
top-down scenarios are highly model-dependent and observational constraints on them are 



-4- 



little at present. Thus we investigate in this study whether the present observations on 
UHECRs can be explained by bottom- up scenarios or not. 

In the bottom-up scenarios, charged particles are accelerated due to Fermi acceleration 
mechanism and obey a power law spectrum ( [Fermi 1949| ). However, the large distances 



between the Earth and potential UHECRs sources like GRBs and AGNs lead to another 
problem called GZK effect ( [Greisen 1966| ; [Zatsepin and Kuz'min 1966|) . As for the protons, 
the energy at which the GZK cutoff takes place (~ 7 x 10^^ eV) is given by the threshold for 
photopion production in the collisions of protons and CMB photons. It is reported that the 
loss length (/loss = Edl/dE) falls below 13 Mpc above 100 EeV ( [Stanev et al. 200q ). As for 
nuclei, the situation is considered to be worse due to the photo-disintegration mechanism 
( [Stanev et al. 2000| ). Thus we assume in this study that the composition of UHECRs is 
proton. 

Taking the GZK sphere for protons into consideration, we can easily understand the 
difficulty of the situation, that is, no plausible source counterparts within the GZK sphere 
have been found within a few degrees from the arrival directions of UHECRs, which are 
considered to be the typical deflection angles of UHECRs (e.g., [Blasi and Olinto 1998| ). 



Moreover, there is a puzzling problem on the distribution of the arrival directions of 



UHECRs (e.g., Blasi and Olinto 1998 ). It is reported that there is no significant large-scale 
anisotropy of arrival direction distribution of UHECRs. On the other hand, there are one 
triplet and three doublets within a separation angle of 2.5 degrees for the 47 cosmic rays 
above 4 x 10^^ eV, and the probability of observing these clusters by a chance coincidence 
under an isotropic distribution is smaller than 1% ( [lakeda et al. 1999| ). We have to give a 
natural explanation for these puzzling problems. 

In this paper, we propose a bottom-up scenario in which the sources of UHECRs 
are correlated with the Super-Galactic Plane (SGP). This is because the probability of 
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observing the clusters by a chance coincidence under such a distribution becomes relatively 
large (~ 10 %) ( [Lemoine et al. 1999|) . In fact, we show in this study that such clusters 



emerge very frequently in our model. It is also shown that the observed energy spectrum 
is well reproduced in our models without any fine-tuned parameter. Moreover, we show 
statistically that the correlation between the arrival directions of UHECRs and SGP can not 
be determined significantly at the present number of data (25 enevts), which is consistent 
with the AGASA data ( [I'akeda et al. 19"^ ). We conclude that about 1000 events, which are 
attained in the very near future as mentioned above, are necessary to determine whether 
the source distribution is correlated with SGP or not. Moreover, we introduce a relatively 
large amplitude of the extra-galactic magnetic field (~ 10 nG) in order to solve the no 
counterpart problem. This is our picture for solving the problems on UHECRs. 

In this study, images of arrival directions of UHECRs are simulated. Optical Redshift 
Survey (ORS) data ([Santiago et al. 1995| ) is used as the source distribution of UHECRs. 
Such a realistic data has not been used in the simulation of propagation of UHECRs in 
the previous works. Inhomogeneity of the arrival directions of UHECRs are also discussed 
statistically by introducing the correlation value, which has not been introduced in the 
previous works (e.g. Lemoine et al. 1999| ). These are new points of this work. We also show 



one realization of arrival directions of UHECRs above 4x10^^ eV in order to compare our 
results with the AGASA data ( [lakeda et al. 1999| ). First and second harmonics analysis 
for our model is also presented in order to compare directly our results with the AGASA 
data ( [Takeda et al. 19991) . Our conclusion has been already stated above. 



In section ^ we show our method of calculation. Results are shown in section 
Summary and discussion are presented in section H. 



2. METHOD OF CALCULATION 
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2.1. Method of Calculation for Propagation of UHECRs 



In this subsection, we describe the method of Monte Carlo simulations for the 
propagation of UHECRs. At first, we assume that the composition of UHECRs is proton. 
We also assume that the initial energy spectrum of the UHECRs obeys the power-law, that 
is, dN/dE oc where N denotes the number of UHECRs. Initial energies of UHECRs 
are assumed to be in the range of (lO^^'^ - 10^^)eV. The number of injected UHECRs in a 
simulation is 10^ except for the cases of Ic = IMpc, where Ic is the correlation length of the 
extra-magnetic field and is explained below. In the case of {B, 1^) = (lOnG, IMpc), the 
number of the injected particles is 10^, where B is the strength of the extra-magnetic field 
and is explained below. In the case of {B, Ic) = (lOOnG, IMpc), the number of injected 
particles is 2.5 x 10^. 

As for the energy loss processes, electron-positron pair creation and photopion 
production in the CMB field are included. Particles below ~ 10^^'^ eV lose their energies 
mainly by pair creations and above it lose their energies mainly by photopion production 
( |Yoshida and Teshima 1993|) . We adopt the formulation of the energy loss rate for the 



pair production on isotropic photons, which has been shown by Chododowski et al. (1992) 
According to them, the energy loss rate of a relativistic nucleus for the pair production on 
isotropic photons is given by 

^4^„icz'!^rd.„(^]^, (1) 



dt rriA J 2 \27/ 

where 7 is the Lorenz factor of the particle, k, = 2k^ ( k is the momentum of the particle 
in units of rriec), n{K), a, tq, Z, and ttla are the photon distribution in the momentum 
space, the fine-structure constant, the classical electron radius, the charge of the particle, 
and the rest mass of the particle, respectively. For the energy range E > 10^^'° eV, can 
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be represented well as 

3 

(p^ nJ^dM (2) 

do ~ -86.07, di ~ 50.96, ^2 ^ -14.45, da = 8/3. (3) 

We also adopt the formulation of the energy loss rate for the photopion production on 
isotropic photons, which has been shown by Achterberg et al. (1999). According to them, 
the interaction length of the photopion production, lp^{Ep), in the CMB can be written as 

where Ep is the proton energy in the cosmic rest frame and £'b is defined as 

Here £th = ?7i^(l + rriT^/'^rnp) ~ 145MeV is the threshold energy of photon in the proton rest 
frame. As for the mean inelasticity, it can be written as 

_ E^ + ml-ml 

= (6) 

where Et is defined as 



Et = + 2mp£o- (7) 

Here Sq is the photon energy in the proton rest frame. As for the maximum spread of 
inelasticity from the mean one, it can be written as 

_ ^{El-m\){E^-ml) 
= 2^ 



= ^{Kp + K^){Kp-K^), (9) 



where m± = nip ± itIt^ and K± = mT^/{mp ^ rriT^). Since photopion productions obey 
the Poisson statistics, we can calculate the energy losses of UHECRs by the photopion 
production using the formulations mentioned above. 
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Next, the details on the extra- magnetic field, which is little known theoretically and 
observationally, is described. Its strength and detailed structure are not known, although 
the strength of the magnetic field in the local super cluster is thought to be relatively 
strong. Only its observational upper limit, ~ IfiG in the meaning of r.m.s., has been 
obtained from the Faraday rotation of the distant sources (|Ryu and Biermann 19981) . In 



this study, 10 nG and lOOnG are adopted as the r.m.s. of the magnetic field. They seems 
to be strong to be sure ( |Sigl et al. 19981 ). However, such strong magnetic field is required 
to solve the problem that there seems to be no plausible source counterparts within the 
GZK sphere within a few degrees from the arrival directions of UHECRs. We also assume 
that the magnetic field is represented as the Gaussian random field with zero mean and a 
power-law spectrum. Thus, {B'^{k)) can be written as 

{B\k)) oc r« for k < fccut (10) 
= otherwise, (11) 

where fccut = 27r//cut = IQtt/Ic characterizes the numerical cut-off scale that is explained 
below and nn is chosen to be -11/3 so as to represent the Kolmogorov spectrum. As for 
the wave numbers of the magnetic field, 8 discrete ones are introduced in this study. In 
practice, ki = 0, i/co, ±2A;o, ±3/co, 4/co are used as wave numbers. Here ko is 2tt/Ic. As for 
the correlation length of the magnetic field, Ic-, three different values are adopted in this 
study. They are 1 Mpc, 10 Mpc, and 40 Mpc. 1 Mpc is widely used as the typical value of 
the correlation length of the extra-galactic magnetic field. 10 Mpc and 40 Mpc represent 
the scale hight and the scale length of the local super cluster, respectively. We separate 
one cubic cell of the size Ic into 512 (= 8^) smaller cells of the size /c/8, which corresponds 
to Zcut mentioned above. The magnetic field are assigned to each small cell. As for the 
boundary condition, the periodic boundary condition is adopted in order to reduce storage 
data for magnetic field components. 
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Finally, we explain the source distribution of UHECRs. In this study, we assume that 
the distribution of sources of UHECRs is proportional to that of the galaxies within the 
GZK sphere. In practice, we use the realistic data from the Optical Redshift Survey (ORS; 
Santiago et al. 1995|) . It is noted that the source distribution will be such like that when 
GRBs occur very frequently in every galaxies. It is also noted that the source distribution 
may be such like that in the top-down scenarios, too. This is because the sources of 
UHECRs such as massive exotic particles will be trapped in the gravitational potential of 
galaxies and/or clusters of galaxies. As for the dependence of our conclusion presented in 
this study on the source distribution, we will report it in the forthcoming paper. 

2.2. Statistics on the Arrival Directions of UHECRs 

2.2.1. Method of Calculation of Correlation Value 



In this study, we introduce the correlation value in order to investigate statistically the 
similarity between the distribution of arrival directions of UHECRs and that of galaxies 
within the GZK sphere. The correlation value, S", between two distributions fi and fs, is 
defined as 

mjs)^ I /^^';^;^ , (12) 

^PUiJilPUsJs) 

where 

Here subscripts j and k discriminate each cell of the sky, Afi(j, k) denotes the solid angle of 
the (j, k) cell, and / means the average of /. In this study, the size of the cell is chosen to 
be 1° X 1°. The meaning of S is as follows. By definition, S ranges from —1 to +1. When 
S = +1, two distributions are same exactly. When S = -1, two distributions are exactly 
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opposite. When S = 0, we can not find any resemblance between two distributions. 

Strictly speaking, the angular images obtained by numerical simulations do not mean 
the exact distributions of the UHECRs which we will detect indeed, but the probability 
density distributions for the arrival directions of UHECRs. Thus, in order to investigate 
the correlation between the results of numerical simulations and the distribution of galaxies 
in the ORS catalogue statistically, we have to make a lot of data samples for the arrival 
directions of UHECRs assuming the numerically obtained probability density distribution. 
In order to estimate the effects of selecting the events stochastically, the number of data 
samples is set to be 1000 for the same number of events of UHECRs, A^, under the same 
condition (e.g., B and Ic)- Hence, a mean value and a standard deviation of Ss can be 
obtained for each N and condition. 

In practice, we choose the 6 different values for the number (A^) of events of UHECRs. 
These are 25, 10^ 320, 10^ 3200, and 10^ 25 is the present number of events of UHECRs 
( |Virmani et al. 2000| ). Other values are chosen arbitrarily only paying attention to the fact 
that the expected number counts detected by the experiment of the next generation such as 
TA, HiRes, Auger, and EUSO are thought to be the order of 10'^ ( [Bhattacharjee and Sigl 
19981). 



2.2.2. Analysis of First and Second Harmonics 



In order to search for the global anisotropy in the arrival directions of UHECRs, we 
apply harmonics analysis to the galactic longitude distribution of events ( [Hayashida et al. 



|1999|) . It should be noted that we do not investigate the right ascension distribution but 
investigate the galactic longitude distribution so as not to suffer from the problem of the 
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incompleteness of the ORS catalog, which does not contain any data within |6| < 20° (see 
Figure ^). We will explain the definition of it below (see Hayashida et al. 1999 in detail). 
The m-th harmonic amplitude, r, and phase of maximum, 9, are obtained for a sample of n 
measurements of phase, 0i, 02, ■ ■ ■, 0n (0 < 0i < 27r) from: 

r = (a2 + 62)1/2 ^^4) 

e = tan'^(6/a) (15) 
where, a = f S"^^ cosm0j, b = sinm0j. 

The following k represents the statistical significance. If events with total number n 
are uniformly distributed in galactic longitude, the chance probability of observing the 
amplitude > r is given by, 

P = exp(-fc), (16) 

where 

k = nr^/4:. (17) 

We take n to be 47 in the range E > 4 x lO^^eV in order to compare our results with 
AG AS A data ([I'akeda et al. l9^ ), although the right ascension distribution of events are 



investigated in their harmonics analysis. We also take n to be 100 in each energy bin in the 
range 10^^-^^ < E < 10^^'^eV. However, it is noted that the results of harmonic analysis for 
the energy range E < 10^^'^eV may not be reliable. This is because the initial energy of 
UHECRs is restricted to be in the range of 10^^-^ < E < lO^^eV to save the CPU time. 



3. RESULTS 



3.1. Arrival Directions and Energy Spectra of UHECRs 
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In this subsection, we show the images of the distribution of the arrival directions of 
UHECRs. We note again that the resolution of the image, that is, the size of the cell is taken 
to be 1° X 1°. It is also noted that these images show the probability density distributions 
of arrival directions of UHECRs. The observed distribution should be interpreted as one 
realization. That is why we have to discuss the correlation among them statistically. In the 
following subsections, the statistical discussion is presented. 

At first, in the left panel of Figure |l], the distribution of the galaxies within 40 Mpc in 
ORS data catalogue is shown. In the right panel of Figure 0, the image in the case of -B = 10 
nG, Ic = 10 Mpc, Ho = 75 km/s/Mpc, and E > 10^^'^ eV is shown. The inter-contour 
interval is 0.5 in the logarithm to base 10 of the integral flux per solid angle. We can find 
that the image is distorted due to the deflection caused by extra-galactic magnetic field. 

Next, the parameter dependence of these images is investigated. In the left panel of 
Figure |^, the image in the same case with Figure |^b but for E > 10^'^ ° eV is shown. We 
find that the image in Figure ^ is more sharp and similar to the image of the distribution 
of galaxies than that in Figure |l]b. It is easily understood because more energetic charged 
CRs propagates more straightly. In the right panel of Figure ^ the image in the same 
case with Figure |]a but for B = 100 nG is shown. Even for the energetic UHECRs, the 
deflection angle becomes relatively large and the image is much distorted as long as the 
amplitude of the extra-galactic magnetic field is taken to be large (~ 100 nG). 

The dependence of these images on the correlation length, l^, is shown in Figure |^. In 
the left panel of Figure ^ the image in the same case with Figure ^ja but for /c = 1 Mpc is 
shown. On the other hand, in the right panel of Figure the image in the same case with 
Figure 0a but for Ic = 40 Mpc is shown. As is shown in Figure ^ the effect of Ic on the 
deflection is relatively weak. This is because the deflection angle is proportional to B and 
/y^, respectively. 



EDITOR: PLACE FIGURE |T| HERE. 
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EDITOR: PLACE FIGURE | HERE. 

In Figure ^, we show one realization of arrival directions of UHECRs above 4x10^^ eV 
in equatorial coordinates for the case of B = lOnG and Ic = IMpc. For comparison, we 
show the distribution of the galaxies within 40 Mpc in ORS data catalogue in equatorial 
coordinates in Figure ^ In Figure ^, the number of events is chosen to be 47 in order to 
compare our results with the AG AS A data ( [lakeda et al. 1999| ). The arrival directions 



of UHECRs are restricted in the range —15° < b < 80° in order to compare our results 
with the AGASA data. Clusters which mean the small-scale anisotropy of UHECRs can 
be found clearly. This result suggests that the probability of observing these clusters is 
relatively high in our model. This is consistent with the conclusion presented by Lemoine 
et al. (1999). In fact, we show in Figure ^ the distribution of the neighbor event (sr~^) as a 
function of the separation angle between the two events whose energies are above 4 x 10^^ 
eV. It can be written as 

where (p denotes the separation angle of the two events. 66 is taken to be 5° in this analysis. 
The number of events in a data sample is chosen to be 47. The number of data samples 
is set to be 1000 in order to obtain the mean value in every bin. We can find clearly the 
significant peak at samll separation angles, which indicates the small-scale anisotropy of 
UHECRs. We also show in Figure the event number of clusters such as doublets and 



- 14- 

triplets within a separation angle of 2.5° in a data sample which contains 47 events, whose 
energies are above 4 x 10^^ eV. The number of data samples is set to be 20 in order to 
obtain the mean values. We can find that the average number of clusters such as doublets 
and triplets is consistent with the AGASA data (P?akeda et al. 1999| ). 



EDITOR: PLACE FIGURE | HERE. 
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Finally, the calculated energy spectra of UHECRs are shown in Figure ||. The 
parameters are set to be -B = 10 nC and Hq = 75 km/s/Mpc. The cases of Ic = 1, 10 and 
40 Mpc are shown respectively. The energy spectra are normalized to the AGASA data 
( [I'akeda et al. l^U^ ) at E' = 100 EeV. We think that the observed energy spectrum is well 
reproduced by all the cases in Figure || without any fine-tuned parameter. 



EDITOR: PLACE FIGURE | HERE. 



3.2. Correlation Value 
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In this subsection, we discuss statistically the correlation between the distribution of 
galaxies within GZK sphere and arrival directions of UHECRs. 

At first, correlation values defined in subsection |2.2.1| are shown in Figure |^ for the 
case of i? = 10 nG, Ic = 10 Mpc and Hq = 75.0 km/s/Mpc indicating the influence of 
the selected energy range. In this flgure, n = 0, n = 1, and n = 2 denote E > 10^^-^eV, 
E > 10^°'°eV, and E > 10^^ '^eV, respectively. It should be emphasized that the present 
observed number (A^) of UHECRs of order 10^ is too small to estimate the flnal correlation 
values. When is close to the order of 10^ - 10^, the correlation values begin to converge 
and flnal correlation values can be estimated. It is well conflrmed in Figure ITUI, in which 



each S is normalized by its flnal value. We emphasize again that a new generation of 
the ground-based large aperture experiments is expected to detect about 1000 UHECRs 
until 2010 (|Zas 200 1|) . Moreover, the mission EUSO is expected to detect about 1000 
UHECRs per year ( phattacharjee and Sigl 1998|) . So we will be able to determine the 



source distribution in the very near future. 



EDITOR: PLACE FIGURE HERE. 
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In Figure [T^, the correlation values are shown in the case of Hq = 75 km/s/Mpc and 
E > 10^° eV indicating the influence of B and Ic- It is conflrmed again that the influence of 
B is greater than Ic, which is indicated in the previous subsection. Furthermore, we can flnd 
again that the number of UHECRs of order 10^ is required to estimate the flnal correlation 
value. 



EDITOR: PLACE FIGURE 
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HERE. 
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3.3. First and Second Harmonics 



We show the results of the first (left panel) and second harmonics (right panel) in the 



galactic longitude in Figure |T2|, |T3|, and |T^. The amplitudes, phases, and chance probability 
of harmonics are shown in these figures. The parameters are set to be i? = lOnG and 
Ic = IMpc. The number of data, n, in one data sample is chosen to be 47 in the energy 



range (4x10^^ - 2x10^° eV) in order to compare our results with the AGASA data (|Takeda 



et al. 19991) . The number of data samples is set to be 1000 in order to obtain the mean 



value and the standard deviation in every bin. 

Compared with the AGASA data, the significant anisotropy on the arrival directions 
of UHECRs can be seen. This may originate from the fact that the results of harmonic 
analysis for the energy range E < 10^^'^eV are not correct. This is because the initial 
energy of UHECRs is restricted to be in the range of 10^^'^ ^ E < lO^^eV to save the CPU 
time. This may also originate from the incompleteness of the ORS database. That is, this 
database does not contain any data within \b\ < 20° (see Figure Also, this problem may 
be solved if the source distribution is slightly changed. For example, this problem may be 
solved if we assume that UHECRs come from some of the galaxies such as AGNs and radio 
galaxies. We will investigate in detail the dependence of our conclusion presented in this 
study on the source distribution in the forthcoming paper. 



4. SUMMARY AND DISCUSSION 



We have performed the detailed numerical simulations on the propagation of the UHE 
protons in the energy range E = (10^^'^ — 10^^'°) eV in the relatively strong extra-galactic 
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magnetic field with strength B = (10, 100) nG within about 40 Mpc. In this case, the 
deflection angles of UHECRs become so large that the no counterparts problem is simply 
solved. 

We have made the images of the angular distribution of UHECRs assuming that the 
source distribution is proportional to the density distribution of matter, or the galaxy 
distribution of the ORS data calatog within 40 Mpc. It is noted that the images obtained 
by numerical simulations do not mean the exact distributions of the UHECRs which we 
will detect indeed, but the probability density distributions for the arrival directions of 
UHECRs. We have found that the influence of Ic on these images are very weak, while that 
of B is very strong. 

We have found many clusters, which mean the small-scale anisotropy, in our simulations. 
This is the advantage to assume that the source distribution is not isotropic. It has been 
also shown that the observed energy spectrum is well reproduced in our models without 
any flne-tuned parameter. 

We have used the correlation value S in order to investigate statistically the similarity 
between the distribution of arrival directions of UHECRs and that of galaxies within the 
GZK sphere. We have found that the values of the parameters are indistinguishable at the 
level of the present number of the events, 25. Also, we have found that we are not still able 
to determine whether the source distribution of UHECRs is correlated with the SGP or not. 

When the number of the detected events becomes O(IO^), the correlation values for 
each parameter set begins to separate and converge. Thus, when the experimental data 
accumulate in the near future, the estimation for the values of parameters will be possible 
by the analysis of the correlation values. The expected number counts by the experiment of 
the next generation such as TA, HiRes, Auger, and EUSO are thought to be the order of 10'^ 
( [Bhattacharjee and Sigl 19981 ). Thus, we will be able to determine the source distribution 
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and values of the parameters in this study in the very near future. 



Compared with the AG AS A data ([lakeda et al. 1999| ), the significant anisotropy on 



the arrival directions of UHECRs can be seen in the analysis of first and second harmonics. 
This may originate from the fact that the results of harmonic analysis for the energy range 
E < 10^^-^eV are not correct. This is because the initial energy of UHECRs is restricted 
to be in the range of 10^^'^ ^ E < lO^^eV to save the CPU time. This may also originate 
from the incompleteness of the ORS database. That is, this database does not contain 
any data within \b\ < 20°. Also, this problem may be solved if the source distribution is 
slightly changed. For example, this problem may be solved if we assume that UHECRs 
come from some of the galaxies such as AGNs and radio galaxies. We will investigate in 
detail the dependence of our conclusion presented in this study on the source distribution 
in the forthcoming paper. 
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Fig. 1. — (a) left panel: angular distribution of galaxies in Optical Redshift Survey (ORS) 
data within 40Mpc. (b) right panel: distribution of arrival directions of UHECRs in the case 
of 5 = 10 nG, Ic = 10 Mpc, Hq = 75 km/s/Mpc, and E > 10^^-^ eV. Galactic coordinate 
is used and left parts correspond to the north galactic hemisphere. The resolution of the 
image is set to be 1° x 1°. The inter-contour interval is 0.5 in the logarithm to base 10 of 
the integral flux per solid angle. 




Fig. 2. — (a) left panel: Same with Figure lb, but for E > 10^°'° eV. (b) right panel: Same 
with Figure 2a, but for B = 100 nG. 
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Fig. 3. — (a) left panel: Same with Figure 2a, but for l^. = 1 Mpc (b) right panel: Same 
with Figure 2a, but for l^. = 40 Mpc 
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Fig. 4. — Arrival directions of UHECRs above 4x10^^ eV in equatorial coordinates for 
the case oi B — lOnG and Ic — IMpc. The number of events is chosen to be 47 in order 
to compare our results with the AGASA data (Hayashida 1999). The arrival directions of 
UHECRs arc restricted in the range —15° < b < 80° in order to compare our results with 
the AGASA data. Clusters which mean the small-scale anisotropy of UHECRs can be found 
clearly. 
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Fig. 5. — Angular distribution of galaxies in ORS data within 40Mpc in equatorial 
coordinate. It is noted that this database does not contain any data within |6| < 20°. 
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Fig. 6. — Distribution of the neighbor event (sr~^) as a function of the separation angle 
between the two events whose energies arc above 4 x 10^^ eV. The number of events in a data 
sample is chosen to be 47. In order to estimate the mean value of every bin, the number of 
data samples is set to be 1000 under the same condition (B — lOnG and Ic — IMpc). 
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Fig. 7. — Number of clusters such as doublets and triplets within a separation angle of 
2.5° in a data sample which contains 47 events, whose energies are above 4 x 10^^ eV. The 
number of data samples is set to be 20 in order to obtain the mean values. 



,v. 26 



Lc=1 
Lc=10 
I c=40 



Mpc 
Mpc 
Mpc 



<n 10^^25 
^ 

< 

?a 

< 

> 

LU 
■□ 




LU 



10^^24 



1 0^23 ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' '~ 

10 100 1000 10000 

Energy [EeV] 

Fig. 8. — Calculated energy spectra of UHECRs. The parameters are set to be 5 = 10 nG 
and Hq = 75 km/s/Mpc. The cases of Ic = 1, 10 and 40 Mpc are shown respectively. The 
energy spectra are normalized to the AGASA data (Hayashida et al. 1999) at E = 100 EeV. 
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Fig. 9. — Correlation values for the case of -B = 10 nG, Ic = 10 Mpc and Hq = 75.0 
km/s/Mpc indicating the influence of the selected energy range, n = 0, n = 1, and n = 2 
denote E > 10^^-^eV, E > lO^^-^eV, and E > lO^^-OeV, respectively. It is noted that the 
present observed number (A^) of UHECRs of order 10^ is too small to estimate the final 
correlation values. When is close to the order of 10^ - 10^, the correlation values begin 
to converge. 
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Fig. 10. — Normalized S by its final value. The adopted parameters are same with Figure 
9. 
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Fig. 11. — Correlation values in the case of Hq = 75 km/s/Mpc and E > 10^° eV indicating 
the influence of B and Ic- It is confirmed that the influence of B is greater than 1^. 



M 
O 



M 
O 



I 



19 



19.5 20 

Log(Energy[eV]) 



20.5 



19 



19.5 20 

Log(Energy[eV]) 



20.5 



Fig. 12. — (a) left panel: amplitude of the first harmonics in galactic longitude in each energy 
bin for the case oi B — lOnG and Ic — IMpc. N is chosen to be 47 in the energy range 
(4x10^^ - 2x10^° eV) in order to compare our results with the AGASA data (Hayashida 
1999). The number of data samples is set to be 1000 in order to obtain the mean value 
and the standard deviation in every bin. (b) right panel: same with Figure 12a, but for the 
second harmonics. 
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Fig. 13. — (a) left panel: same with Figure 12a, but for the phase of the first harmonics, 
(b) right panel: same with Figure 13a, but for the second harmonics. 
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Fig. 14. — (a) left panel: same with Figure 12a, but for the chance probability of the first 
harmonics, (b) right panel: same with Figure 14a, but for the second harmonics. 



